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Abstract 

For  portable  fuel  cell  systems  a  multitude  of  applications  have  been  presented  over  the  past  few  years.  Most  of  these  applications  were 
developed  for  indoor  use,  and  not  optimised  for  outdoor  conditions.  The  key  problem  concerning  this  case  is  the  cold  start  ability  of  the  polymer 
electrolyte  membrane  fuel  cell  (PEMFC).  This  topic  was  first  investigated  by  the  automotive  industry,  which  has  the  same  requirements  for 
alternative  traction  systems  as  for  conventional  combustion  engines. 

The  technical  challenge  is  the  fact  that  produced  water  freezes  to  ice  after  shut-down  of  the  PEMFC  and  during  start-up  when  the  temperature 
is  below  0  °C. 

To  investigate  the  basic  cold  start  behaviour  isothermal,  potentiostatic  single  cell  experiments  were  performed  and  the  results  are  presented. 

The  cold  start  behaviour  is  evaluated  using  the  calculated  cumulated  charge  transfer  through  the  membrane  which  directly  corresponds 
with  the  amount  of  produced  water  in  the  PEMFC.  The  charge  transfer  curves  were  mathematically  fitted  to  obtain  only  three  parameters 
describing  the  cold  start-up  with  the  cumulated  charge  transfer  density  and  the  results  are  analysed  using  the  statistical  software  Cornerstone 
4.0. 

The  results  of  the  statistic  regression  analyses  are  used  to  establish  a  statistic -based  prediction  model  of  the  cold  start  behaviour  which 
describes  the  behaviour  of  the  current  density  during  the  experiment.  The  regression  shows  that  the  initial  start  current  mainly  depends  on  the 
membrane  humidity  and  the  operation  voltage.  After  the  membrane  humidity  has  reached  its  maximum,  the  current  density  drops  down  to 
zero.  The  current  decay  also  depends  on  the  constant  gas  flows  of  the  reactant  gases. 

Ionic  conductivity  of  the  membrane  and  charge  transfer  resistance  were  investigated  by  a  series  of  ac  impedance  spectra  during  potentiostatic 
operation  of  the  single  cell  at  freezing  temperatures.  Cyclic  voltammetry  and  polarisation  curves  between  cold  start  experiments  show 
degradation  effects  by  ice  formation  in  the  porous  structures  which  lead  to  significant  performance  loss. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Sub  zero  temperature  operation  of  portable  fuel  cell 
systems  is  a  relatively  new  topic  and  not  much  literature  is 
published.  Cargnelli  et  al.  operated  a  50  W  PEMFC  system 
in  a  climatic  chamber  at  —  42  °C  and  used  for  the  start-up 
a  catalytic  burner  to  heat  up  and  humidify  the  reaction  air 
[1].  Chu  et  al.  operated  a  50  W  fuel  cell  system  at  —  10°C 
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for  9h  [2].  Datta  and  Velayutham  reports  of  a  500  W  fuel 
cell  system  which  was  operated  without  any  problems  in  the 
Antarctic  [3]. 

More  fundamental  in  situ  investigations  of  non- 
isothermal,  galvanostatic  single  fuel  cell  experiments  below 
0  °C  were  performed  by  Kagami  et  al.  At  constant  gas  flow  it 
is  described  that  the  cathode  electrode  surface  is  reduced  by 
ice  formation  and  therefore  the  reaction  is  inhibited.  With  the 
empirical  assumption  that  the  cathode  surface  is  diminished 
proportional  with  the  amount  of  frozen  water,  simulations 
were  done  which  describe  the  voltage  decay  over  time  in 
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Nomenclature 

BET 

Brunauer-Emmett-Teller  porosity  analysis 

c 

exponent  variable,  constant 

CV 

cyclic  voltammetry 

EDS 

energy  dispersive  X-ray  spectroscopy 

F 

Faraday  constant  (96485  C  mol-1) 

GDL 

gas  diffusion  layer 

i 

current  density  (mA  cm-2) 

linear, stari 

current  density  at  start  of  slow  decrease 
(mA  cm-2) 

Jpeak 

peak  current  density  during  cold  start 
(mA  cm-2) 

Ltart 

current  density  at  initial  cold  start  (mA  cm-2) 

MEA 

membrane  electrode  assembly 

mn2o,prod  area  specific  product  water  (g  cm  2) 

Mw2o 

molar  weight  of  water  (g  mol-1) 

MPL 

microporous  layer  between  electrode/GDL 

P 

influence  parameter  in  regression  polynomial 

P450 

reference  power  density  at  450  mV 
(mW  cm-2) 

R 2 

stability  index 

*adj 

adjusted  stability  index 

Ret 

charge  transfer  resistance  from  Nyquist  plot 
(m<2  cm2) 

^m,c 

membrane/contact  resistance  from  Nyquist 
plot  (m£2  cm2) 

SEM 

scanning  electron  microscopy 

^q,cum(0 

cumulated  charge  transfer  density  (C  cm-2) 

fitted  cumulated  charge  transfer  density 
(mC  cm-2) 

^q,max(0 

max.  cumulated  charge  transfer  density 
(C  cm-2) 

Sq, peak, end  cum.  charge  transfer  density  at  end  of  peak 
(mC  cm-2) 

5q, peak, start  cum.  charge  transfer  density  at  start  of  peak 
(mC  cm-2) 

^63 

time  when  Sq^t(t)  reaches  63%  of  Sq,mSLX  (min) 

^air 

air  flow  rate  (ml  min-1) 

V 

voltage  (mV) 

Voc 

open  circuit  voltage  (mV) 

X[ 

value  of  influence  parameter 

y 

response  variable  in  regression  polynomial 

Zstart 

Impedance  at  1000  Hz  before  cold  start-up 
(mQ) 

good  correlation  with  experimental  results.  Beneficial  for  the 
cold  starts  are  low  current  densities  which  prevent  voltage 
decay  by  freezing  of  product  water  [4,5]. 

In  this  paper  we  present  the  results  of  potentiostatic  single 
cell  cold  start-up  measurements  under  isothermal  conditions 
at  — 10  °C.  The  current  density  is  increasing  very  fast  after  ini¬ 
tial  start-up  and  later  it  decays  towards  zero.  This  behaviour 
is  ascribed  to  freezing  of  product  water  in  the  cathode.  To 


find  out  the  main  influencing  operational  parameters  on  the 
cold  start  behaviour,  statistic  methods  are  used.  Finally,  a 
statistic-based  prediction  model  of  the  current  decay  during 
isothermal  cold  start  was  developed. 

To  get  a  deeper  insight  into  the  dynamic  processes  in 
the  cathode  during  cold  start-up  dynamic  electrochemical 
impedance  spectroscopic  analyses  were  done,  which  show 
that  the  membrane/contact  resistance  as  well  as  the  charge 
transfer  resistance  is  changing  with  the  amount  of  produced 
water. 

Impedance  spectroscopy  at  80  °C  was  also  used  by  Cho  et 
al.  to  find  out  changes  of  single  fuel  cell  characteristics  by  ex 
situ  thermal  cycling.  The  author  concludes  that  the  contact 
resistance  in  the  fuel  cell  is  increasing  after  thermal  cycles 
because  of  worse  contact  between  the  membrane  and  the 
electrode,  whereas  the  membrane  ionic  conductivity  itself 
is  not  affected  [6] .  The  impact  of  ex  situ  thermal  cycling  on 
the  degradation  of  Nation®  membranes  was  examined  by 
McDonald  et  al.  It  is  reported  that  the  ionic  conductivity  and 
the  mechanical  properties  are  hardly  changing.  However, 
opening  up  of  the  molecular  structure  of  the  polymer  in 
connection  with  enlargement  of  hydrophilic  areas  was  found 
[7].  Cappadonia  et  al.  investigated  Nation®  membranes  by 
thermal  cycling  and  it  was  found  that  two  different  water 
environments  exist  below  0°C.  Water  phase  transition 
depends  on  the  water  content  in  the  membrane.  Higher  water 
contents  lead  to  larger  pore  diameters  in  the  membrane 
and  therefore  to  a  lower  freezing  temperature  of  water. 
In  water-rich  samples  phase  transitions  were  observed  at 
260  K  [8]. 

Nevertheless,  performance  degradation  can  be  observed, 
as  described  in  [6,9].  A  degradation  rate  of  about  2.8%  at 
80  °C  and  0.6  V  per-freeze-thaw  cycle  down  to  — 10  °C  is 
observed  as  well  as  a  reduction  of  about  6%  per  thermal  cy¬ 
cle  of  the  electrochemical  active  surface  area  of  the  cathode. 
It  was  found  that  the  pore  size  distribution  in  the  electrode 
changes  to  larger  pores  by  thermal  cycles. 

Our  measurements  with  cylic  voltammetry  (CV)  also 
show  degradation  of  the  cathode  electrode  surface  by  isother¬ 
mal  sub  zero  operation,  and  furthermore  changes  in  hy- 
drophobicity  of  both,  the  microporous  layer  (MPL)  and  the 
gas  diffusion  layer  (GDL)  on  the  cathode  side.  We  finally 
also  found  performance  degradation  at  450  mV  and  30  °C  of 
more  than  5%  per  each  cold  start  experiment. 


2.  Experimental  and  results 

2.7.  Maximum  cumulated  charge  transfer  density 

The  cold  start  experiments  were  carried  out  with  a  single 
cell  with  an  active  area  of  46  cm2 .  The  flowfield  had  a  double 
serpentine  structure  and  the  cell  assembly  was  done  using 
commercial  standard  components  consisting  of  GDL,  MPL 
and  catalyst-coated  membrane. 
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Before  freezing  the  fuel  cell  a  reference  performance  mea¬ 
surement  at  450  mV  with  dry  hydrogen  (stoichiometry  1.5) 
and  dry  air  (stoichiometry  3.0)  at  30  °C  was  recorded.  This 
reference  power  density  is  called  £>450.  Each  reference  point 
was  measured  at  least  20  min  whereby  gas  flows  were  ad¬ 
justed  automatically  to  the  actual  current. 

After  this  reference  operation  the  gases  were  switched  to 
dry  nitrogen  (N2),  and  the  anodic  as  well  as  the  cathodic  side 
of  the  fuel  cell  were  purged  with  1.01  min-1  at  30  °C.  The 
purging  process  continued  until  the  high  frequency  resistance 
at  1000  Hz  increased  by  a  factor  of  about  100  and  more.  Then, 
N2  flow  was  switched  off  and  the  cell  was  cooled  down. 

After  drying  and  freezing  of  the  cell,  the  start-up  experi¬ 
ments  were  carried  out  potentiostatically.  Gas  flow  rates  as 
well  as  the  cell  temperature  were  kept  constant.  The  dew  point 
temperatures  of  the  gases  are  between  —  10°C  and  —  5°C. 
When  the  graphite  plates  reach  the  required  temperature,  dry 
hydrogen  and  dry  air  are  switched  on.  After  open  circuit  volt¬ 
age  has  arisen,  potentio static  load  is  applied  and  the  current 
over  time  is  recorded.  For  evaluation  of  the  cold  start  experi¬ 
ments,  the  cumulated  charge  transfer  density  SqjCum(0  serves, 
which  is  calculated  by  (1). 

5q,cum(0  —  f  1  dt  (1) 

JO 

5q,cum(0  integrates  over  the  dynamic  current  density  and  it  di¬ 
rectly  correlates  with  the  amount  of  product  water  raH2o,prod, 
given  by  Eq.  (2). 


m  h2o  = 


’q,cum 


(0 


IF 


^h2o 


To  find  out  the  main  influencing  operating  parameters  15 
experiments  were  carried  out  at  constant  temperature  of 
— 10  °C.  In  Table  1,  all  varied  operational  parameters  are 
listed. 

After  the  fuel  cell  is  dried  out,  cooled  down  to  — 10  °C  and 
Voc  has  arisen,  potentiostatic  load  is  applied.  The  resulting 
current  reacts  over  time  always  with  the  same  principle  be¬ 
haviour,  that  is  a  fast  increase  which  reaches  within  the  first 
2  min  its  maximum  and  then  decays  towards  zero.  To  pre¬ 
dict  the  maximum  cumulated  charge  transfer  density  Sq?max, 
a  mathematical  fitting  function  was  found  which  parame- 
terises  5q,cum(0-  This  exponential  function  is  calculated  by 

(3). 

Sq.lit(0  =  5q-max(l  -  (3) 


The  quality  of  a  cold  start  experiment  is  evaluated  by 
5q,max-  A  linear  regression  analysis  was  performed  to  pre¬ 
dict  the  main  influencing  operational  parameters  on  Sq> max 
using  the  statistic  software  Cornerstone®  (Brooks  Automa¬ 
tion).  It  fits  a  response  y,  e.g.  S^max,  with  a  polynomial  term 
consisting  of  a  constant  c;  the  values  jq-  of  the  p  influence 
parameters,  e.g.  applied  voltage  V,  multiplied  by  constant 
factors  at.  In  (4),  a  general  polynomial  equation  is  shown 
where  linear  effects  are  taken  into  account  (at  xf),  interaction 
effects  between  each  parameter  (ay  Xy)  and  quadratic  effects 

(an  Xj 2). 

p  p  p 

y  =  c  +  °ixi  +  E  a‘jXii  +  E  anF  (4) 

i~  1  i,j=l,i^j  i= 1 

The  fit  quality  of  the  regression  polynomial  is  evaluated  by 
the  stability  index  R 2  and  the  adjusted  stability  index  R2d ■. 

R 2  is  the  relation  between  the  statistical  spread  of  the  model 
and  the  statistical  spread  of  the  real  measured  values.  For  test 
series  with  a  low  number  of  experiments  R2d ■  is  better  suited, 
because  it  accounts  additionally  for  the  number  of  experi¬ 
ments  and  the  number  of  influence  parameters.  The  values 
of  R2  and  R2d ■  are  between  0  and  1 .  For  significant  regres¬ 
sion  model  the  stability  index  should  be  near  one  [13].  Using 
the  calculated  polynomial  which  represents  the  regression 
model,  it  is  possible  to  predict  5q,max  within  the  regression 
boundaries. 

Fig.  1  illustrates  the  regression  result  for  Sq,max  as  pre¬ 
dicted  response  graph.  Three  parameters  were  found  which 
have  a  significant  impact  on  the  value  of  S^max.  The  largest 
effect  has  the  reference  power  density  77450,  which  depends  on 
the  degree  of  degradation.  The  influence  of  the  air  flow  rate 
ua u-  and  the  membrane  humidity  indicated  by  Zstart  is  lower. 
The  coefficient  correlation  between  these  three  parameters  is 
lower  than  0.4,  which  indicates  that  they  are  independent.  On 
the  y-axis  the  predicted  response  from  the  linear  regression 
polynomial  is  shown  as  the  function  of  the  values  of  £>450 
(here  adjusted  to  250mWcm-2),  u air  (600  ml  min-1)  and 
Zstart  (40  mQ).  The  predicted  value  for  Sq,max  is  47.99  C  cm-2 
with  a  standard  deviation  of  ±12.9  C  cm-2. 

This  regression  result  indicates  that  two  operational  pa¬ 
rameters  can  be  adjusted  by  the  operator  to  achieve  a  high 
charge  transfer  below  0°C.  The  membrane  should  be  dry, 
indicated  by  a  high  Zstart  and  the  air  flow  rate  should  be  also 
as  high  as  possible.  The  third  and  most  strongly  influencing 
parameter  is  the  overall  performance  parameter  77450.  With 


Table  1 


Varied  operating  parameters  for  start-up  experiments  at  —  10°C 


Operating  parameter 

Symbol/unit 

Values 

Applied  voltage  for  cold  start-up 

y/mV 

200/400/600 

Air  flow  rate 

w  air /ml  min-1 

147-1000 

Hydrogen  flow  rate 

WHo/mlmin-1 

20-158 

Cell  impedance  after  nitrogen  purging 

Zpurging/mf2 

70-1500 

Cell  impedance  after  applying  gases,  directly  before  cold  start 

Zstart/mf2 

7-50 

Power  density  at  450  mV,  30  °C,  dry  gases,  stoichiomety  1. 5/3.0 

P45o/mW  cm-2 

210-350 

410 
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250  600  40 


p_450  u_air  Z_start 

Fig.  1.  Predicted  Sq,max  as  function  of  the  dominating  three  operational  pa¬ 
rameters.  Linear  regression  of  all  operational  parameters  (/?2dj  =  0.87). 

higher  performance  the  local  electrode  temperature  is  higher 
[10]  and  it  is  assumed  that  crystallisation  of  water  starts  in 
the  more  distant  areas  of  the  GDL  and  not  immediately  in  the 
electrode.  Therefore,  a  longer  operation  is  possible  before  ice 
formation  hinders  the  reaction. 

2.2.  Current  decay 

To  analyse  the  initial  starting  behaviour  at  — 10  °C  the 
current  density  i  is  referred  to  the  cumulated  charge  transfer 
density  Sq,Cum(0-  For  statistical  evaluation  it  is  necessary  to 
parameterise  these  current  density  curves,  which  is  done  by 
taking  five  characteristic  curve  points  into  account. 

To  find  out  which  operational  parameters  affect  the  shape 
of  the  current  decay  curve,  linear  regression  analyses  were 
made  for  each  point.  Interactions  of  the  main  parameters  (V, 
^air,  mh2?  Zstart,  P45o)  are  included.  The  goodness  of  the  re¬ 
gression  analysis  is  indicated  by  R2d j  which  is  for  all  points 
between  0.855  and  0.998,  except  for  point  5  OS^max),  where 
R2d j  is  0.640.  Table  2  shows  all  coefficients  for  the  statistic 
prediction  polynomial. 

As  an  example  the  regression  model  calculates  Sqjmear 
(Point  4)  with  a  polynomial  expression  given  by  the  constant 
c  and  coefficients  multiplied  by  the  adjusted  values  of  the 


Fig.  2.  Statistical  based  fit  of  current  decay  curves  at  —10  °C  and  different 
operating  conditions  as  shown  in  legend.  Experimental  data  are  indicated 
with  open  symbols  (A,  0,  □,  O)  and  the  belonging  statistical  5-point  fit 
with  filled  symbols  (A,  ♦,  ■,  •). 

influence  parameter,  in  this  example  u air  and  Zstart: 

Sq, linear  =  3.9557  -  0.0024^air  +  0.4184Zstart 
T 0 . 0004 W  air  Z  start 

In  Fig.  2,  four  experimental  results  for  the  current  density 
over  the  charge  transfer  density  for  different  operating  condi¬ 
tions  are  shown.  Additionally,  the  fitting  curves  of  the  above- 
introduced  statistical  5 -point  fitting  model  are  included.  This 
demonstrates  that  it  is  possible  to  describe  the  current  den¬ 
sity  decay  over  the  cumulated  charge  transfer  density  by  the 
operational  parameters. 

The  statistic-based  model  estimates  the  dominating  factors 
for  each  curve  point.  In  Table  3,  the  percentage  influence  is 
shown. 

The  regression  shows  that  the  initial  start  current  /start 
mainly  depends  on  the  membrane  humidity  Zstar t  and  the 


Table  2 


Coefficients  for  polynomial  prediction  model 


Point 

© 

(D/(D 

© 

© 

Ltart 

5q,  peak,  start 

*peak 

5q,  peak,  end 

5q,  linear 

^linear 

5q,max 

(mA  cm-2) 

(mA  cm-2) 

(mA  cm-2) 

(mA  cm-2) 

(C  cm-2) 

(mA  cm-2) 

(mA  cm-2) 

c 

-977.8981 

-22.4874 

144.9494 

7.4862 

3.9557 

-1.9906 

-26.0127 

V 

2.9383 

-0.2209 

-0.0180 

-0.0056 

^air 

0.9819 

0.0229 

-0.4419 

0.0104 

-0.0024 

0.0582 

0.0247 

MHt 

0.6911 

-0.6934 

Zstart 

-36.3194 

0.3027 

-3.1301 

0.3954 

0.4184 

0.4772 

P450 

3.8487 

0.1402 

0.1094 

0.1663 

^air^m 

Vhair 

-0.0014 

-0.0008 

-0.0001 

FZstart 

0.0158 

0.0036 

P45()V 

Zstart 

-0.0104 

-0.0078 

0.0107 

^air-Zstart 

-0.0123 

0.0004 

0.0004 

P450  wair 

P450  x  Zstart 
t/WHo 

0.0888 

0.0021 

0.0009 
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Table  3 


Estimated  influence  of  main  parameters  on  current  decay  curves 


Point 

© 

l  start 

©A3) 

© 

© 

5q,max 

Sq,  peak,  start 

*peak 

^q, peak, end 

Sq,  linear,  start 

flinear, start 

V 

34% 

— 

40% 

25% 

— 

59% 

— 

^air 

19% 

68% 

92% 

5% 

13% 

20% 

5% 

— 

69% 

— 

4% 

— 

— 

— 

Zstart 

43% 

36% 

32% 

87% 

85% 

— 

26% 

P450 

12% 

— 

34% 

— 

— 

10% 

26% 

operation  voltage  V.  Humidification  of  the  membrane  and 
the  anode  electrode  sets  in  by  product  water  until  the  current 
density  reaches  its  peak  plateau  /peak-  The  membrane  reaches 
maximum  protonic  conductivity  at  Sq,peak,end-  Now  product 
water  starts  to  decrease  the  surface  of  the  porous  electrodes 
and  the  current  density  drops  down  very  fast.  The  level  to 
which  it  drops  down  /linear, start  is  strongly  dominated  by  the 
operating  voltage  V.  Lower  voltage  is  connected  with  higher 
waste  heat  production  at  the  electrode  resulting  in  a  higher 
/linear, start  because  product  water  freezes  in  larger  distance 
from  the  electrode. 

This  interpretation  is  supported  by  the  fact  that  all  charge 
transfer-related  statements  are  highly  influenced  by  the  dry¬ 
ing  status  of  the  cell  whereas  all  current-related  statements 
can  be  interpreted  through  heat  production  and  are  signifi¬ 
cantly  influenced  by  the  cell  performance  ^450,  the  cell  volt¬ 
age  V  and  the  air  flow  u air. 

Blockages  of  the  gas  channels  by  accumulated  product 
water  droplets  would  lead  to  higher  pressure  drops  over  the 
flowfield  at  constant  gas  flows.  This  is  verified  by  compar¬ 
isons  of  the  pressure  drop  during  fuel  cell  operation  with  and 
without  (at  Voc)  current  production. 

It  is  assumed  that  during  cold  start-up  experiments  with 
constant  gas  flows  and  stable  pressure  drop,  most  of  the  prod¬ 
uct  water  does  not  reach  the  flowfield  channels,  but  freezes 
in  the  porous  structures  of  the  electrode,  MPL  and  the  GDL, 
which  corresponds  with  the  proposed  mechanism  presented 
in  [4,5]. 


Fig.  3  shows  an  isothermal  cold  start  experiment  at 
— 10  °C,  constant  gas  flow  rates  of  150  ml  min-1  air  and 
120  ml  min-1  hydrogen.  Ten  EIS  were  applied  in  potentio- 
static  mode  at  600  mV  with  an  effective  excitation  of  ±5  mV. 
The  current  density  increases  during  the  first  three  EIS  within 
1.5  min  from  16  up  to  38  mA cm-2.  EIS  No.  4  and  5  were 
taken  during  the  fast  current  decay  between  2.5  and  7.5  min 
after  starting  the  experiment.  The  measuring  time  for  the  fol¬ 
lowing  EIS  No.  6  to  EIS  No.  10  was  nearly  80  min,  which 
means  that  the  current  decayed  much  slower  and  that  it  was 
more  stable  during  EIS  recording. 

Fig.  4  shows  10  Nyquist  plots  of  measured  EIS  during 
this  cold  start  experiment.  Next  to  the  last  measuring  point 
the  applied  lowest  frequency  is  printed.  Assuming  the  simple 
electrode  model  of  the  cathode  as  suggested  in  [11]  consist¬ 
ing  of  the  membrane/contact  resistance  Rm  c  in  series  to  the 
parallel  charge  transfer  resistance  7?ct,  and  the  double-layer 
capacity  Cd,  one  can  estimate  the  values  for  Rm  c,  which  is  ad¬ 
dressed  to  the  left  basis  point  of  the  semicircle  and  Rct ,  which 
represents  the  diameter  of  the  semicircle  in  the  Nyquist  plot. 
The  semicircles  of  EIS  No.  1  and  No.  5  show  the  above- 
mentioned  artefact  in  the  last  measuring  point  because  the 
cell  current  was  too  unstable  while  EIS  was  recorded. 

By  linear  regression  analyses  the  most  important  influence 
parameters  on  the  decrease  of  Rm  c  at  the  initial  start  of  the  ex¬ 
periment  (5q?Cum  <  10  C  cm-2)  were  identified.  In  Fig.  5,  the 
result  shows  that  Rm  c  decreases  significantly  with  increasing 


2.3.  Electrochemical  impedance  spectroscopy 

A  single  fuel  cell  with  a  double  serpentine  flowfield  and  an 
active  area  of  33  cm2  was  used  for  electrochemical  impedance 
spectra  (EIS)  during  isothermal,  potentiostatic  cold  starts. 
Different  frequency  ranges  were  measured,  applying  an  ef¬ 
fective  voltage  excitation  of  ±5  mV  (sine  wave).  The  anode 
side  (H2)  served  as  the  reference  electrode  and  was  connected 
to  the  counter  electrode.  The  cathode  side  (air)  was  connected 
to  the  working  electrode.  This  set-up  is  also  described  in  [6]. 
The  upper  frequency  was  1000  Hz  and  the  lower  frequency 
varied  between  0. 1  and  0.04  Hz,  because  of  the  high  dynamic 
behaviour  of  the  current  density  over  time.  Since  the  current 
decay  is  fast,  impedance  spectra  must  be  measured  in  a  very 
short  time  otherwise  the  results  show  artefacts  (fold  back  of 
the  semicircle)  at  low  frequencies,  as  the  measuring  time  is 
inversely  proportional  to  the  applied  frequency. 
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Fig.  3.  Current  density  decay  and  operating  voltage  (600  mV  ±5  mV)  of 
isothermal  cold  start-up  at  — 10  °C,  150  ml  min-1  air  and  120  ml  min-1 
H2.  The  voltage  signal  is  modulated  by  the  applied  EIS  frequency 
(1000-0.1  Hz/0.04  Hz)  with  an  effective  excitation  of  ±5  mV. 
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Fig.  4.  Series  of  10  EIS  during  a  cold  start-up  experiment  at  600  mV, 
lSOmlmin^1  air  and  120  ml  min-1  H2  flow.  Frequency  range  of  EIS  be¬ 
tween  1000  Hz  and  40  mHz  (lowest  frequency  is  printed  next  to  the  last 
point).  The  voltage  amplitude  is  ±5  mV. 


120  3  400 


Fig.  5.  Predicted  membrane/contact  resistance  Rm^c  (m£2)  as  function  of  the 
dominating  three  parameters.  Linear  regression  of  all  experimental  parame¬ 
ters  (R2ad-  =  0.75). 

Sq,cum-  Higher  operating  voltage  and  hydrogen  flow  lead  to 
lower 

While  there  is  a  dependency  between  Rm  c  and  the  hy¬ 
drogen  flow  one  can  find  a  significant  influence  of  the  air 
flow  for  the  charge  transfer  resistance  Rct ,  which  is  mainly 
dominated  by  Sq.  Fig.  6  illustrates  the  regression  result  for 
Rct  in  a  predicted  response  graph  for  all  Sq,Cum  higher  than 
10C  cm"2. 

Statistical  evaluation  of  AC  impedance  spectra  during 
isothermal  cold  start-up  shows  a  strong  dependency  be¬ 
tween  the  cumulated  charge  transfer  density  Sq,Cum,  the 
membrane/contact  resistance  and  the  charge  transfer 
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Fig.  6.  Predicted  charge  transfer  resistance  Rc{  (m£2)  as  function  of  the  dom¬ 
inating  three  operational  parameters.  Linear  regression  of  all  operational 
parameters  ( R^  =0.83). 


Fig.  7.  Charge  transfer  resistance  Rct  and  membrane/contact  resistance  Rm,c 
as  function  of  the  cumulated  charge  transfer  density. 

resistance  Rct.  Fig.  7  illustrates  the  typical  behaviour  of  both 
resistances  during  a  cold  start  experiment. 

Two  effects  support  the  shape  of  the  current  decay  curves. 
At  the  initial  start  Rm  c  is  high  and  stabilises  later  at  a  much 
lower  level.  As  a  result,  the  ohmic  losses  improve.  This  re¬ 
sult  confirms  the  formerly  formulated  assumption  that  the 
membrane  gets  humidified  by  product  water  even  at  sub  zero 
temperatures.  Rct  reacts  converse;  with  longer  experiment  du¬ 
ration  it  increases  and  therefore  the  current  density  decays. 
This  effect  is  addressed  to  a  reduction  of  the  active  surface 
of  the  cathode  electrode  by  ice  formation. 

2.4.  Degradation  effects 

During  several  series  of  cold  start  experiments  strong 
performance  degradation  at  450  mV  was  found.  By  cyclic 
voltammetry  on  the  cathode  side  a  decrease  of  the  hydrogen 
desorption  peak  after  10  cold  start-up  experiments  was 
observed  as  shown  in  Fig.  8.  The  anode  side  (flushed  with 
H2)  served  as  the  reference  electrode  and  was  connected  to 
the  counter  electrode.  The  cathode  side  (flushed  with  N2) 
was  connected  to  the  working  electrode.  After  Vqc  dropped 


Fig.  8.  Cyclic  voltammograms  after  cold  start-up  experiment  1  and  10  with 
a  scanning  rate  of  50  mV  s- 1 . 
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Fig.  9.  Degradation  of  performance  72450,  surface  area  Aec  and  increase  of 
membrane/contact  resistance  Rmc  after  10  cold  start  experiments. 


below  100  mV  and  stayed  stable,  the  cell  was  cycled  three 
times  between  0  and  1.4  V  with  a  scan  rate  of  50  mV  s_1 . 

The  electrochemical  active  area  Aec  can  be  estimated  us¬ 
ing  the  area  of  the  hydrogen  desorption  peak,  the  scan  rate 
and  the  specific  charge  for  a  hydrogen  monolayer  on  plat¬ 
inum,  which  is  210  jxC  cm-2.  The  fractional  coverage  of  the 
hydrogen  monolayer  of  0.77  as  proposed  by  Breiter  et  al.  [12] 
is  not  included. 

Fig.  9  illustrates  all  detected  degradation  effects  during 
10  cold  start  experiments  for  Aec,  the  reference  power  den¬ 
sity  P450  and  the  high  frequency  membrane/contact  resistance 
Rm  c  at  the  beginning  of  the  cyclic  voltammogram. 

The  decrease  in  the  calculated  electrochemical  active  sur¬ 
face  Aec  from  about  50  to  40  m2  g-1  is  connected  with  higher 
activation  polarisation  [6]  and  the  increase  in  the  high  fre¬ 
quency  resistance  from  230  to  270  mQ  cm2  does  explain 
higher  ohmic  losses. 

The  average  percentage  changes  per  cold  start-up  are  cal¬ 
culated  in  Table  4.  The  values  are  in  the  same  range  as  re¬ 
ported  in  [6,9]  and  show  that  it  is  absolutely  necessary  to 
improve  cold  start  and  shut-down  procedures  to  prevent  per¬ 
formance  degradation. 


Table  4 


Degradation  by  cold  start  experiments 


Value 

Unit 

Change  per  cold  start  experiment 

P450 

mW  cm-2 

-5.4% 

Ac 

2  -1 
m  g 

-2.4% 

An,c 

m£2  cm2 

+1.9% 

After  disassembling  the  cell  it  was  found  that  the  flow- 
field  channel  pattern  is  visible  on  the  cathode  GDL  as  dark 
lines  (microscope  image  in  Fig.  10a,  magnification  3x).  A 
hydrophobicity  test  showed  that  the  channel  pattern  is  much 
more  hydrophilic  (Fig.  10b)  and  water  droplets  agglomerate. 

Ice  formation  could  deform  the  pore  structure  in  the  MPL 
and  the  GDL.  Oh  et  al.  found  by  BET  analysis  after  ther¬ 
mal  cycles  that  the  pore  size  distribution  of  the  electrode 
changed  to  larger  pores  [6,9].  The  same  effect  could  occur  in 
the  MPL  and  the  GDL.  Other  effects  could  also  play  a  role, 
like  washed-out,  hydrophilic  polymer  chains  of  the  mem¬ 
brane  or  loss  of  PTFE  in  the  GDL. 

To  prove  the  loss  of  hydrophobicity  in  the  GDL  we 
performed  contact  angle  measurements  with  10  pul  water 
droplets.  Fig.  11  shows  that  the  contact  angle  on  the  dark 
areas  is  about  112°,  whereas  on  the  not-coloured  areas  the 
contact  angle  is  about  133°.  On  the  anode  side  the  same  hy¬ 
drophobicity  test  shows  a  uniform  distribution  and  no  channel 
pattern.  The  anode  contact  angle  is  136°  and  comparable  with 
hydrophobic  parts  of  the  cathode  GDL.  This  indicates  again 
that  mainly  the  cathode  is  affected  by  ice  formation  during 
cold  starts. 

In  this  cell  assembly  also  microporous  layers  were  used, 
which  are  embedded  by  the  MEA  and  the  GDL.  After  dis¬ 
assembly  of  the  cell  it  was  found  that  the  hydrophobicity  of 
the  cathode  MPL  has  only  changed.  As  well  as  on  the  cath¬ 
ode  GDL  water  accumulated  on  the  channel  structure  after 
the  spray  test.  In  Fig.  12,  the  channel  pattern  on  the  cathode 
MPL  is  visible  by  the  agglomeration  of  water  droplets.  The 
image  was  taken  under  a  microscope  with  a  magnification  of 
about  3  x . 


Fig.  10.  Flowfield  channel  pattern  on  cathode  GDL  after  cell  disassembly  and  spray  test  with  water.  The  hydrophilic  channel  pattern  is  visible. 
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Fig.  11.  Contact  angle  measurements  with  10  pi  water  droplets,  (a)  uncoloured  areas  of  cathode  GDL  (0=  133°),  (b)  coloured  channel  area  of  cathode  GDL 
(0  =112°)  and  (c)  anode  GDL  (0  =  136°). 


Fig.  12.  Flowheld  pattern  on  the  microporous  layer  of  the  cathode.  Water 
droplets  cumulate  in  the  more  hydrophilic  channel  structure. 

Further  investigations  to  clarify  the  structural  changes 
will  be  done  by  mercury  porosimetry,  BET  analysis  and 
SEM/EDS  analyses. 


3.  Conclusions 

By  isothermal  potentiostatic  cold  start  measurements  of 
single  fuel  cells  under  different  operating  conditions  we 
found  that  water  freezes  in  the  porous  structures  of  the  cath¬ 
ode  electrode,  microporous  layer,  and  the  GDL.  Statistic  eval¬ 
uation  of  the  experiments  showed  that  dryer  membranes  and 
high  air  flow  rates  are  beneficial  for  high  charge  transfer  rates. 
Furthermore,  it  was  found  that  performance  degradation 
plays  an  important  role  on  the  charge  transfer  too.  As  pre  stage 
to  physical  models  we  presented  statistic-based  models  which 
explain  the  characteristics  of  the  current  density  during  cold 
start-up  as  function  of  all  operating  parameters  and  estimates 
their  significant  influence.  From  these  results  we  conclude 
that  product  water  increases  first  the  membrane  humidity 
even  at  — 10  °C,  and  after  it  has  reached  its  maximum,  product 
water  floods  the  porous  structures  of  the  electrode,  MPL  and 
GDL,  and  freezes.  Therefore,  the  current  density  decays.  This 


assumption  is  proved  by  electrochemical  impedance  spec¬ 
tra  during  isothermal  potentiostatic  cold  start  experiments, 
which  show  that  the  membrane/contact  resistance  decreases, 
whereas  the  charge  transfer  resistance  increases  over  time. 
For  both  values  the  charge  transfer  and  the  amount  of  prod¬ 
uct  water,  respectively,  plays  a  significant  role.  Ice  formation 
leads  to  a  strong  degradation  of  the  cell  performance  and 
the  electrochemical  active  surface  area  of  the  cathode.  The 
membrane/contact  resistance  increases  with  the  number  of 
cold  start-up  experiments  and  structural  changes  on  the  cath¬ 
ode  GDL  and  MPL  were  observed,  leading  to  changes  in 
hydrophobicity.  To  prevent  degradation,  elaborated  start-up 
and  shut-down  procedures  have  to  be  developed  in  the  future 
to  overcome  the  obstacle  cold  start-up  of  PEM  fuel  cells. 
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